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The propagation of s- and p-polarized light through quasiperiodic multilayers, consist- 
ing of layers with different refractive indices, is studied by the transfer matrix method. 
In particular, we focus on the transmission coefficient of the systems in dependency on 
the incidence angle and on the ratio of the refractive indices. We obtain additional bands 
with almost complete transmission in the quasiperiodic systems at frequencies in the 
range of the photonic band gap of a system with a periodic alignment of the two ma- 
terials for both types of light polarization. With increasing incidence angle these bands 
bend towards higher frequencies, where the curvature of the transmission bands in the 
quasiperiodic stack depends on the metallic mean of the construction rule. Additionally, 
in the quasiperiodic systems for p-polarized light the bands show almost complete trans- 
mission near the Brewster's angle in contrast to the results for s-polarized light. Further, 
we present results for the influence of the refractive indices at the midgap frequency of 
the periodic stack, where the quasiperiodicity was found to be most effective. 
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1 Introduction 

The insight in the photonic properties of materials is important for the construction of cus- 
tomized new optical devices. In this context, quasicrystals, which are regarded to have a 
degree of order intermediate between crystals and disordered systems, are of special interest 
because they often possess unexpected physical properties due to their complex symmetries 
[H El IS]- This qualifies them for the potential application in several optical devices such 
as single-mode light-emitting diodes, polarization switches, and optical filters [UG3E1- An 
interesting application is also the construction of microelectronic devices that are based on 
photons rather than on electrons, so that they potentially can be the electromagnetic analogue 
to semiconductors [8] . 

The theoretical study of photonic properties of one-dimensional systems is based on the 
transfer matrix method and the concept of aperiodic mathematical sequences, as e.g. the 
Fibonacci sequence, the Thue-Morse sequence, or Cantor sequences [HI [9l [TUl E] • Such one- 
dimensional systems can be relatively easily produced in reality and a comparison of the 
theoretical and the experimental results shows a good agreement [T2| [T3]. 
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Figure 1: Transmission of light through a quasiperiodic stack consisting of the materials A 
and B with different refractive indices. The light wave is incident on the surface 
under an angle Parts of the light wave are transmitted through the stack (E^) 
and other parts are reflected (E®). 



While in literature most results are only given for zero incidence angle, s-polarized light 
and fixed refractive indices, we compare the transmission of s- and p-polarized light through 
different multilayers in dependency on the incidence angle of the light and the refractive 
indices. 

In particular, the configuration of the layers of the so-called metallic-mean quasicrystals 
can be constructed by the inflation rule 

[b ^ A 

v ' = \a^a-b • (1) 

iterated m times starting from the symbol B. We refer to the resulting sequence after m 
iterations as the mth-order approximant with the length f m given by the recursive equation 
fm — fm-2 + afm-i with the initial values fo = f\ = 1. Depending on the parameter a the 
inflation rule generates different metallic means, i.e., the lengths of two successive sequences 
satisfy the relation 

lim -p^- = lim 4? = T ( a ) ' ( 2 ) 



fm-1 m ^°° f] 



where r(a) is an irrational number with the continued- fraction representation r(a) = [a] = 
[a, a, a, ...]. Also, the number of layers A and number of layers B in an approximant 
are related by this metallic mean. 

For example, a = 1 yields the well-known Fibonacci sequence related to the golden mean 
[1] = (l + y/E)/2, the case a = 2 corresponds to the octonacci sequence with the silver 
mean tap- = [2] = 1 + y/2, and for a = 3 one obtains the bronze mean tb z = [3] = (3 + \/l3)/2. 
In general the relation r(a) = (a + \/a 2 + 4)/2 holds [T^ [TT]. 

The outline of this paper is as follows: In Sec. [2] we give an introduction to the transfer 
matrix method used for the calculations and in Sec. [3] we present our results and discuss them. 
This is followed by a brief conclusion in Sec. [H 
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2 Transfer Matrix Method 



The propagation of light through a layered system as shown in Fig. Q] is commonly investi- 
gated by transfer matrix methods [9l [TOj [15l [16] . Regarding the geometry we compare the 
propagation of linearly polarized light with the electric field perpendicular to the plane of 
the light path (s-polarization) and with the electric field vector lying in the plane of the light 
path (p-polarization) (cp. Fig. ED [T7] . 

Within one layer the electric field E = Ee y can be described as a superposition of a right 
and a left traveling plane wave E = E J r e lk i x ~ lujt + ^ e -ikjx-iut ^ w } iere fa. — n .fc denotes the 
wave number of the light in the jth layer with refractive index nj, k the wave number in the 
vacuum, and uj the frequency of the light. 




Figure 2: Electric and magnetic field vectors E and H at the interface between layer j and 
j + 1 with different refractive indices n a and np for s-polarized light (left) and 
p-polarized light (right). 



The propagation of the waves through the stack is given on the one hand by their refraction 
at the interfaces between the layers and on the other hand by their propagation through the 
layers. Addressing the first point, the boundary conditions at the interfaces between the layers 
require the tangential component of the electric field E and the magnetic field H = x E to 
be continuous across the boundaries [IB]. Applying the (£"+, £'_) T notation used by Kohmoto 
et al. P [TO] with E+ = E T + E\ and E- = —i (E T — E\) , the interface matrices T a p = 
for s- and p-polarized light, respectively, 

T a(3 = ( n np cosflfl ) = | np cos^ j (3) 

\ riccos^c/ \ n a cos ^ J 

describe the transfer of the light wave from medium a to medium /?. Thereby, the incidence 
angle the emergence angle tip, and the corresponding refractive indices n a and np are 
related by Snell's law according to n^sin^Q, = npsiw&p. Further, total reflection occurs at 
incidence angles $ > i?totai — arcsin (n^/n^). 

On the other hand, the propagation of the light waves within the layers results in a phase 
difference, which is comprised in the propagation matrix 

T = ( COS ^ sm ^7^j u\ 
7 siny? 7 cos(/? 7 y 
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The phase difference for the wave length A = 2tt/ k of the light and a layer thickness d 7 amounts 
to (fj = 2kd 1 n 1 cos # 7 [T7]. For the choice of the wave length of the incident light Ao = 2tvc/(jJq, 
we use the commonly applied quarter wave length condition tla^a — ^b^b — Ao/4, which 
leads for an incidence angle i? = to an identical optical wave path of the light in the two 
materials A and B of the stack [T5J El E] • Hence, we obtain 

(fij — ^-7rcos^ 7 — — 7rcosi? 7 . (5) 

The transfer matrix of the overall system relates the incident light E®, the reflected light 
E®, and the transmitted light E^ (cp. Fig. []]) by 

E+ \ , ^Qcl ( Ea 



E -),r \ e -)« <6) 

and is given as a combination of the different interface and propagation matrices according 
to a certain quasiperiodic sequence. In general, the recursive equation 

mZ = {mZ-,TmZ. 2 (7) 

is applicable for m > 2 with Mq a = TabTbTba and Mf a = T4. For refractive indices of the 
surrounding medium rtc 7^ tla the corresponding transfer matrices Tqa and T^c have to be 
added at the edges. 

In particular, we are interested in the calculation of the transmission coefficient T (also 
known as transmittance) of the light through the stack defined as 

^_ |^| 2 _ 4(detM) 2 



|£?| 2 |M| 2 + 2detM ' 

where |M| 2 corresponds to the sum of the squares of all four matrix elements of M. Further, 
for the complete stack it can be shown that det M = 1. 



3 Results 

In this section we comprise the results for several metallic-mean quasicrystals. Thereby, we 
focus on the change of the transmission coefficient T in dependence on the inflation rule V, 
the light polarization, the ratio of the refractive indices u — ua/tib, and the incidence angle 
1). Here we only show results for one system size f(C a ). In general, increasing the number 
of layers results in a reduction of the width of the transmission bands and the occurrence 
of new transmission bands in a selfsimilar manner [HI HH] • In the limit of infinite systems 
the transmission bands form a Cantor set with Lebesgue measure 0, and it is assumed that 
either only complete transmission or complete reflectance occurs [21 [TO]. However, for physical 
applications usually systems with relatively small numbers of layers are used 0,11]. 

In Fig. E] the results for the transmission coefficients T(#, uS) are shown in dependence on 
the reduced frequency oj/ojo and the incidence angle $ for s-polarized and p-polarized light. 
For 7? = the transmission spectrum is periodic with respect to the frequency. However, 
for all angles 1? > the photonic transmission bands bend towards higher frequencies for 
increasing incidence angles and this effect becomes stronger for higher frequencies. The main 
difference of the two light polarizations occurs in the region around the Brewster's angle 
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Figure 3: Transmittance T(i?, uj) of light through a stack of the materials A and B, which are 
arranged according to the octonacci sequence C^ g with /§ = 577 layers. Results 
are shown in dependence on the reduced frequency uj/ujo versus the incidence angle 
i? for via — Tic, u — 0.5 and (a) s-polarized and (b) p-polarized light. The blue 
and red lines indicate the expected bending of the transmission bands according to 
Eqs. (1TTJ) and (fT3]) . The green lines show the regions of constructive interference 
according to Eq. (j9j). 



$Br = arctantx , which corresponds to an incidence that produces an angle of 90° between 
the reflected and refracted ray. For p-polarized light we obtain in this region much higher 
transmission coefficients and the different transmission bands emerge more clearly even for 
large incidence angle, where we hardly see any transmission in the s-polarized case. Of course, 
for incidence angles i? = $Br — 1-107 one always obtains a transmission coefficient T = 1 for 
u — 0.5 and p-polarized light. 

In Fig. |4]the transmittance T($, uj) is displayed for different metallic-mean sequences. Com- 
paring with the results for a periodic alignment of the layers A and £>, one can see that there 
is a typical photonic band gap at uj/ujq = 6+1/2 (be No) for the periodic chain for # = 0. 
For the quasiperiodic sequences this band gap appears to be much wider, but is intersected 
by an increasing number of new lines of moderately high transmission for increasing a of the 
inflation rule. For uj/ujq = b the characteristics of the transmittance T(??,a;) hardly change 
for the different construction rules. The behavior at the photonic band gap is in consistency 
with the observation that for uj/ujq — 6+ 1/2 = 0) the quasiperiodicity is most effective [9]. 
Again we observe the occurrence of almost complete transmission coefficients for p-polarized 
light near the Brewster's angle ??Br — 0.464. 

Comparing Figs. [3]and[4j i.e., the transmission for u > 1 and u < 1, we observe for u > 1 
a relatively large range of frequencies with nearly complete transmission for large incidence 
angles and for u < 1 certain frequencies with almost complete transmission over a large 
range of angles. These regions correspond to constructive interference for a layer B which 
is embedded in medium A, i.e., the transmission coefficient of the matrix TabTbTba yields 
1 for these parameter setups. This is satisfied for (fs = (n E Z) and yields a bending 
according to 

uj n n 

— . (9) 

ujo cosi?^ cosarcsin(^ _1 sin??) 

Note that for ua — tlq the relation <&a — # holds. Hence, these regions of total transmission 
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Figure 4: Transmittance T(#, u) of different quasiperiodic stacks for s-polarized (upper row) 
and p-polarized light (lower row) and a ratio of the refractive indices u = 2: (a),(e) 
a periodic sequence (AB) n with / = 500 layers, (b),(f) the Fibonacci sequence Cy± 
with fu = 610 layers, (c),(g) the octonacci sequence C$ g with f$ = 577 layers and 
(d),(h) the bronze-mean quasiperiodic sequence Cf T with f$ = 469 layers. The red 
and green lines have the same meaning as in Fig. [3l 



in Figs. [3] and H] bend independently of the underlying (quasiperiodic or periodic) alignment 
of the B layers as indicated by the green lines. 

A similar interference argument can be applied to estimate the curvature of the transmission 
bands in the quasiperiodic stack. This behavior originates from the difference of the optical 
light paths in the layers A and B for incidence angles i? > 0. For the derivation we assume 
that the transmission coefficients T(cj,$ = 0) and T(uj' \$ > 0) yield almost the same results 
if the overall phase difference for the stack is equal. At first, we consider only the behavior 
arising from the change of the phase difference within one layer. Hence, for layer B the 
relation (^^(cj,^ = 0) = (fB^^B > 0) has to be fulfilled, which yields a behavior according 
to 

f ^0) _ ^ = 0) 

cos v b cos (arcsm [u sin v) ) 
Likewise, the relation for layer A can be derived with 

i "(* = Q) . 

cosv 
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Figure 5: Limit behavior of the bending of the transmission bands with the incidence angle $ 
for quasiperiodic stacks in dependency on the metallic mean r(a): the green, blue, 
and red lines correspond to Eqs. (jTU]h (ITT]) , and (p3J) respectively. 



The actual bending of the transmission bands is obtained by a superposition of Eqs. (fTUl) 
and (fTTl) . For the quasiperiodic stacks the overall phase difference depends on the numbers 
of layers A and B in the stack according to Eq. (pQ). This results in the limit of an infinite 
stack in 

f T~~ — ¥A + — ¥B J fm (12) 



1 + T 1 + T 



and yields a bending according to 



r cos v + cos (arcsm [u sin ?/)) 

These functions are shown in Fig. [5] for the cases u — 0.5 and u = 2. We obtain a bending 
between the behaviors arising from one single layer of medium A or B. In particular, the 
bending increases with the parameter r (respectively a) in the case u < 1 and shows the 
reverse behavior for u > 1. Further, for the quasiperiodic systems the bending according to 
Eq. ([TS]) approaches for r — » 1 the bending of a periodic stack and for r ^> 1 the behavior 
arising from one single layer B. 

Comparing Eq. (IT3]) with the numerical results of the transmission for the quasiperiodic 
stacks for different values of r (cp. Fig. [4]), we find a very good resemblance of the bending 
behavior for small and intermediate angles. While this indicates that the curvature of the 
transmission bands is solely caused by the phase difference, which is described by the prop- 
agation matrix, the differences in the transmission intensity for s- and p-polarized light near 
the Brewster's angle can only be caused by the different transfer matrices. 

In the following, the characteristics of the transmittance T are investigated with respect to 
the change of the ratio of the refractive indices u. In Fig. E]we study T(oj,u) in dependency 
on the reduced frequency cj/ujo and the ratio of the refractive indices for the golden and silver 
mean model. Again it is clearly visible that the quasiperiodicity has the largest effect for 
uj w cjo/2. In this region the transmittance T(cj, u) strongly varies with the used construction 
rule. By increasing the parameter a, more and more bands appear and the width of the 
transmission bands becomes smaller. In the periodic case complete transmission occurs for 
this frequency only for the obvious case u = 1, whereas for the quasiperiodic systems there 
are several bands in this region with a high transmission coefficient. Again, the bending of the 
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Figure 6: Transmittance T(uj,u) of quasiperiodic stacks for $ = 0.4 and (a),(b) s-polarized 
and (c),(d) p-polarized light: (a),(c) the Fibonacci sequence Cf 4 u with f u = 610 
layers and (b),(d) the octonacci sequence C% g with fg = 577 layers. The red and 
green lines have the same meaning as in Fig. El 

bands is well described by Eq. (IT3]) . Further, for p-polarized light we obtain almost complete 
transmission near the ratio of the refractive indices u^ T = 1/ tan 7? for $ = 0.4, i.e., in this 
case the Brewster condition is met. This leads to a wide range of values u E [1, u^ r ] for which 
almost complete transmission can be observed. 

The results for the influence of the refractive indices on the incidence angle $ are shown 
in Fig. where the red lines in Fig. [3(c), (d) indicate the Brewster's angle. A comparison of 
the plots for s- and p-polarized light shows that the transmission is again strongly increased 
around the Brewster's angle. Comparing the results for the golden and silver mean model, we 
observe a similar structure of the transmittance T($, u) for small angles and near the regions 
of total reflection. While for small incidence angles transmission in a periodic system is only 
possible for values of u close to 1, in the considered quasiperiodic systems there is a whole 
range of possible ratios of refractive indices u with almost complete transmission. 

4 Conclusion 

We studied the transmission of light through quasiperiodic multilayers in dependency on the 
underlying structure, the incidence angle, and the light polarization. We obtained additional 
bands with almost complete transmission in the quasiperiodic systems at frequencies in the 
range of the photonic band gap of a periodic system for both types of light polarization. We 
found that with increasing incidence angle these bands bend towards higher frequencies and 
derived a relation in dependency on the metallic mean r(a). This indicates that the bending 
is solely caused by the overall phase difference occurring during the passage of the stack. 
Further, in the quasiperiodic systems for p-polarized light the bands showed almost complete 
transmission near the Brewster's angle $Br m contrast to the results for s-polarized light. In 
quasiperiodic systems almost complete transmission occurs for a significant range of ratios u 
in contrast to the periodic case. 
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Figure 7: Transmittance T($, u) for the midgap frequency uj — and the same systems as 
in Fig. El For parameters in the white areas in the bottom right of the plots total 
reflection occurs. In plots (c),(d) the Brewster condition is indicated by a red line. 
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